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With the use of Gutzwiller-projected variational states, we study the renormalization of the 
current carried by the quasiparticles in high-temperature superconductors and of the quasiparticle 
spectral weight. The renormalization coefficients are computed by the variational Monte Carlo 
technique, under the assumption that quasiparticle excitations may be described by Gutzwiller- 
projected BCS quasiparticles. We find that the current renormalization coefficient decreases with 
decreasing doping and tends to zero at zero doping. The quasiparticle spectral weight Z+ for 
adding an electron shows an interesting structure in k space, which corresponds to a depression of 
the occupation number rik just outside the Fermi surface. The perturbative corrections to those 
quantities in the Hubbard model are also discussed. 

PACS numbers: 71.10.-w 



I. INTRODUCTION 

In recent years, it has been acknowledged that 
ground-state properties of high-temperature supercon- 
ductors may be reasonably well described with the 
help of Gutzwiller-projected wave functions How- 
ever the main challenge of any candidate theory of 
high-temperature superconductivity is the description of 
finite-temperature properties such as the superconduct- 
ing transition and the pseudogap phenomenon in under- 
doped cuprates. One of the first issues related to the 
finite-temperature physics of high-temperature supercon- 
ductors is the structure of low-lying excitations. Within 
the framework of Gutzwiller-projected wave functions, 
the first steps in studying the excitations have been re- 
cently made: the quasiparticle spectrum and the quasi- 
particle spectral weight have been calculated 0, Q ■ In 
our paper, we complement the previous studies with the 
analysis of the current carried by the quasiparticles. The 
magnitude of the quasiparticle current has a direct physi- 
cal implication in reducing the superfluid density at finite 
temperature, which eventually determines the supercon- 
ducting transition temperature in the underdoped regime 
0, 0. Furthermore, the deviation of the quasiparticle 
current from the prediction of the BCS theory may pro- 
vide a helpful insight in the physics of high-temperature 
superconductivity. 

The reduction of the superfluid density ns{T)/m by 
thermal quasiparticles at the nodes of a d-wave super- 
conductor has been computed by Lee and Wen as 

n'nsiT) _h^nM 21n2 ^vp.. ^ ... 



where vp and wa are the velocity of the nodal quasipar- 
ticles perpendicular and parallel to the underlying Fermi 
surface, and a is the phenomenological Landau param- 
eter 6] which renormalizes the current carried by the 
quasiparticle 

j(k) = -ea^F. (2) 

Experimentally, ns{T)/m can be related to the London 
penetration depth A, and the ratio vf/va may be ex- 



tracted independently from a thermal-conductivity mea- 
surement . This makes a an experimentally accessible 
quantity for a variety of doping values H, . 

In the first part of the paper, we focus on computing 
the current renormalization a for different doping values. 
We find that it decreases with decreasing doping, and 
that it is roughly constant along the Fermi surface at 
all dopings. The contribution of particles and holes to 
the total quasiparticle current allows us to picture the 
"effective Fermi surface" where the electron contribution 
crosses over to the hole contribution. We observe that 
this "effective Fermi surface" deviates considerably from 
the original Fermi surface of the unprojected BCS state. 
This reveals the particle-hole asymmetry produced by the 
Gutzwiller projection. 

In the second part of the paper, we discuss another 
renormalization parameter: the quasiparticle spectral 
weight Z+ for adding an electron. The momentum de- 
pendence of Z+ shows a pocket structure at the diagonal 
of the Brillouin zone just outside the Fermi surface. We 
further discuss the relations and bounds on Z+ in the t-J 
model, as well as corrections arizing from the rotation to 
the Hubbard model. 

For our analysis, we take the minimal two dimensional 
model for strongly interacting electrons on a lattice, the 
Hubbard model. Following the usual procedure (see, e.g., 
Refs. 0, ITo|) - we first study the wave function for its 
strong-coupling limit, the t-J model and then include the 
first-order correction in t/U due to the doubly-occupied 
sites. The Hamiltonian for the t-J model for our system 
is defined on the two-dimensional square lattice by 

+J^(S, -Sj -n,n,/4), (3) 

with t/J ~ 3. Here is the electron creation operator 
at site i with spin a ~ (|, |). The (i, indicates nearest 
neighbors, = ^4 Q,aa,/3Ci,/3, and Ui = cj„Ci,Q. The 
Hamiltonian lO is then supplemented with the constraint 
that no double occupancy is allowed. 
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We use the variational Monte Carlo technique to calcu- 
late expectation values of operators given our trial wave 
functions. For the t-J model, we consider two re- 

lated trial wave functions: the ground state wave function 
5'gs} and the wave function for the excited state |5'ex)- 
For the ground state, we use the Gutzwiller-projected 
d-wave singlet. 



|*Gs) = PdPn\^bcs) 



(4) 



where Pd — [1 — n^^riij] is the Gutzwiller pro- 
jection operator onto the subspace with no doubly 
occupied states, and Pm is the projection operator 
onto the subspace with N particles. I^bcs) = 

Ilk 1 + OkCk |C^k i 1*^)' where ^ is the Fourier trans- 
form of real-space electron creation operator c\ ^ . Fol- 
lowing the standard BCS definitions for a d-wave singlet 
state, 



Ok 

Ak 



ek + 

— ^var(cOS kx 



k 

cos ky ) 



(5) 

(6) 
(7) 



with ek = — 2(cosA:a; -I- cosfcy). Not only has this d-wave 
Gutzwiller-projected wave function been shown to give 
good variational energies for the t-J model relative to 
other possible phases, but also it has correctly reproduced 
many properties of the superconducting state. 0,01 

For the trial wavefunction of the low-lying excited 
states, we take the natural ansatz of the Gutzwiller- 
projected Bogoliubov quasiparticle 



l^-Exlk, cr)) = PdPniIJ'^bcs) 



(8) 



Since the overall normalization of the wave function is 
of no importance, we may also rewrite the trial excited 
state as 



|*Ex(k,a)) = Pz5P^vcLl*BCs). 



(9) 



Throughout the paper, we suppress the k and a variables 
on I^'ex) for notational convenience. The expectation 
values of any operator O in the variational ground and 
excited states will be often denoted as (0)03 and (C')ex 
respectively. 

In our simulations, we use the optimal values of Avar 
and Uyar calculated for the ground state trial wave func- 
tion [l^, both for the ground state and for the ex- 
cited state. These values minimize the expectation value 
(T^tj) GS of the physical t~J Hamiltonian at a fixed con- 
centration of holes. 

We assume boundary conditions that are anti-periodic 
in the x direction and periodic in the y direction so that 
we avoid the singularity in ak along the nodal diagonal, 
(0,0) to (7r,7r). A drawback of this choice of boundary 
conditions is that we are unable to calculate quantities 
exactly on the nodal diagonal, and instead calculate ex- 
pectation values for nearby k points. This becomes an 
important source of error when we compute expectation 
values at the nodal point. We try to lower this error by 



looking at larger systems in order to get better resolution; 
however, we are limited by our computing resources. 

This trial excited state and various related ones have 
been studied previously by other groups 0, 0, 0, 0| . 
The energy dispersion of the low-energy quasiparti- 
cles has been found to be of the BCS type E(k) = 

1 /2 

[^^ -I- A^] , but with renormalized values of the gap 
and bandwidth. The nodal point obviously coincides 
with the nodal point of the unprojected wave function 
^BCS and slowly shifts inwards from (7r/2,7r/2) along 
the diagonal of the Brillouin zone as the hole doping in- 
creases. 



II. CURRENT CARRIED BY THE 
QUASIPARTICLES 

In this section, we investigate the current carried by 
quasiparticles as a function of their momenta and doping. 
The current carried by the excited state |\l/Ex(k, tr)) is 
defined as 



4 



0) 



EX, 



(10) 



where {ij) represents a sum over all links. Note that the 
orientation and direction of the link determines its con- 
tribution to the vector j. We can interpret the excited 
state trial wave function as the ground state for a system 
of N particles to which we add one unpaired electron so 
that our state has a net spin and charge. The Gutzwiller 
projection enforces strong correlation between the added 
electron and the A^-electron ground state, so that we ex- 
pect an effective quasiparticle with renormalized param- 
eters. 

In the standard BCS theory (without the Gutzwiller 
projection), the current carried by the quasiparticles is 
jk — evk where Vk = dsk/dk is the velocity from the 
underlying normal state and not the velocity from the 
quasiparticle dispersion dE/dii. This is the result of the 
fact that in BCS theory the excitations are a superpo- 
sition of particle and hole states and that these states 
carry opposite charge but also move in opposite direc- 
tions. The underlying metallic state can have some Fermi 
liquid correction to the current carried by the particles 
and holes. This correction is then carried through to 
the quasiparticle current in the superconducting state, 
jk = aevk. In the case of the Gutzwiller projected trial 
excited state that we are studying, we see that current 
does still approximately follow the shape of the disper- 
sion of the underlying metal and that the quasiparticle 
current renormalization a can be calculated from the ra- 
tio of j/vk. 



A. Current as a function of k 

First we examine j, the current carried by the quasi- 
particle, as a function of k for various dopings. At most 
dopings, we find a distribution with a structure similar 
to what one would expect from the tight binding model. 
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FIG. 1: The top middle and bottom figures are respectively 
plots of the current j, j|, and j| as a function of the wavevec- 
tor k for a 12x12 system with 13 holes, x — 0.09. The vectors 
are drawn starting at the k-point at which the current is cal- 
culated, their length is proportional to the current magnitude. 



In the top plot of Fig. a typical example of the current 
carried by the quasiparticle as a function of wavevector 
is plotted. We compare the direction of the current as 
a function of momentum to both the quasiparticle dis- 
persion, dE(]i.)/dk, and to the underlying dispersion of 
the normal state Vk. We find that the shape of the cur- 



ie - (0,0) to (7t,3l) 

FIG. 2: The magnitude of the current |j| measured in units 
of t along the nodal diagonal, (0, 0) to (n, tv). These runs were 
done on the 10x10 system with the doping of 0.01, 0.05, 0.09 
and 0.17 (increasing magnitude). 



rent indeed approximately follows the dispersion of the 
normal state and not the dispersion of the quasiparticles. 
This is the same as in the BCS theory so the Gutzwiller 
projection does not change this aspect of the physics. 

For intermediate doping values, the magnitude of the 
current reaches its maximum value near the center of 
the Brillouin zone, and therefore near the nodal point. 
However, for strongly underdoped simulations, x < 0.05, 
we find that the maximum moves inward along the nodal 
direction, becoming closer to (7r/4, 7r/4) as the doping 
approaches zero. (See Fig. |2J) 

We also look at how much of the total current is be- 
ing carried by the up spins and down spins individually. 
By restricting the J^a equation ^| so that we con- 
sider a and a =i separately, we can investigate this 
property of our trial excited state. Since we are now dis- 
tinguishing between these two spins, it is important to 
note that we define our trial wave function as adding an 
up-spin to the system. Given this, we find that the cur- 
rent of our state is almost entirely carried by either the 
up spins or the down spins depending on whether or not 
we are inside or outside of the effective Fermi surface as 
can be seen in the lower two plots of Fig. Q Inside of the 
Fermi surface, all of the current is carried by the down 
spins and outside of the Fermi surface all of the current 
is carried by the up spins. This is the qualitatively the 
same as in the unprojected case where the up and down 
spin currents have factors of and respectively. 

To make a more detailed comparison to the BCS the- 
ory, we define the "particle contribution" to the current 
as 



jjW ■ jtof(k) 
|jtot(k)|2 



(11) 



where jtot = Jx + Jt is the total current carried by the 
quasiparticle. In the BCS theory, nj(k) = |wkP, it 
takes the values between zero and one, and the isoline 
rij(k) = 0.5 coincides with the Fermi surface. In Fig. |21 
we show the contour plots of nj(k) for different values 
of doping, together with the Fermi surface for the corre- 
sponding unprojected BCS wave functions. We see that 
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FIG. 4: Magnitude of the nodal current, |jnodai| measured 
in units of t, plotted as a function of doping. The solid lines 
are the data from a 10x10 lattice and the dashed lines from 
a 20x20 one. The lower curves are for the trial wavefunc- 
tion for the t-J model and the upper curves are for the trial 
wavefunction for the Hubbard model (to the lowest order in 
t/U = 1/12). 



uated at k = (O.Stt, 0.47r) and for a 20x20 system at 
k = (0.457r, O.Stt). In Fig. ^ we see that there is agree- 
ment to within the error between the data calculated for 
two lattices of different size, and we expect that the ac- 
tual nodal current would also be within these errors. 

In the 20x20 system, we can study the doping as low 
as 0.005 (2 holes), and our results indicate that the cur- 
rent apparently decreases down to zero with decreasing 
doping. 



FIG. 3: The particle contribution to the quasiparticle current 
nj(k) defined in Eq. lllll . The top and the bottom plots 
correspond to 13 and 31 holes in the 12x12 system (doping 
X = 0.09 and x = 0.22, respectively). The contour lines are 
nj(k) = 0.1, 0.3, 0.5, 0.7, and 0.9 (from left to right). The 
thick dashed line is the Fermi surface of the unprojected state, 
the thick solid dot marks the position of the node. Small 
crosses indicate the positions of data points (we have used 
the reflection about the x-y diagonal to double the density of 
data points). 

the "effective Fermi surface" defined by nj(k) does not 
follow the original Fermi surface of the BCS state, but 
bends outwards in the (0, tt) regions. Thus it effectively 
acquires an inward curvature similar to the effect of the 
negative t' hopping term. 



B. Current as a function of doping 

We noted earlier that the magnitude of the current 
has a maximum near the nodal point and that this max- 
imum decreases as a function of doping. In Fig. ^ we 
plot the magnitude of the nodal current jnodai = j(k = 
(7r/2, 7r/2)) versus doping for both 10x10 and 20x20 sys- 
tems. Because of our choice of boundary conditions, we 
do not calculate the current at the true nodal point. 
For a 10x10 system, the "nodal point" is actually eval- 



C. Rotation to the Hubbard model 

So far we have studied the properties of fully projected 
wave functions. We would now like to extend our simu- 
lations from the t-J model back to the Hubbard model. 
This can be done in the standard way by employing a uni- 
tary transformation e*"^ that decouples the Hilbert space 
of the Hubbard Hamiltonian so that there are no matrix 
elements connecting those subspaces with different num- 
bers of doubly occupied sites. Following MacDonald et 
gLjigj. we determine this transformation as a power se- 
ries in {t/U), so that the subspaces are decoupled order 
by order. To the first order, the rotation generator is 
given by 

*5 - ^ [Ti - , (12) 

where Ti and T_i are defined to be the parts of the ki- 
netic energy operator that increase and decrease, respec- 
tively, the number of doubly occupied sites by one. 

With the use of this rotation, we can replace comput- 
ing the expectation value of any operator O in the ground 
state of a Hubbard model by computing the expectation 
value of the rotated operator {e^^'Oe^^^} in the ground 
state of the t-J model. We use the same variational wave 
function for the t-J model as described above and com- 
pute the lowest-order correction (linear in t/U) to the t-J 
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FIG. 5: The renormalization of the quasiparticle current a = 
Ijnodail/i'F as a function of doping. The open circles are for 
runs on a 10x10 lattice and the solid ones are for runs on a 
20x20 one. The solid line is the fit to the 20x20 data by Eq. 
lITCl at X* = 0.09. 



expectation value. This procedure has been applied, for 
example, in the work of Paramekanti et al. 10] for cal- 
culating the expectation value of the occupation-number 
operator. 

We compute the rotation correction to the value of 
the nodal current and find that it does not qualitatively 
change its doping dependence. The corrected value of 
the current is plotted as the upper curves in Fig. 0] (for 
t/U= 1/12) 



D. Quasiparticle current renormalization a 

Beyond just studying the current itself, we are particu- 
larly interested in the quasiparticle current renormaliza- 
tion factor a. As we noted earlier the current and the 
slope of the dispersion are coUinear within the error of 
our simulations, so we can define a = j/vk. 

We look at the quasiparticle current renormalization 
parameter a as a function of doping. We are interested 
in a for the lowest lying excitations, i.e., for those at 
the nodal point. Although we have calculated the Fermi 
velocity at the nodal point, more precise data for this 
velocity as a function of doping are available from the 
work by Yunoki et al. Ref. fl. We use those Fermi ve- 
locity data in conjunction with our nodal current results 
to calculate a at the nodal point. 

In Fig. [51 we plot the nodal value of a as a function of 
doping. We find that a goes to zero at zero doping. 

It is interesting to compare our results with the pre- 
dictions of slave-boson theory. In this theory, the low ly- 
ing excitations are x bosons which carry charge with an 
effective hopping matrix element proportional to t and 
fermions which carry spin with an effecting hopping pro- 
portional to J. These excitations are coupled to gauge 
fluctuations. When the gauge fluctuations are treated at 
the Gaussian level, we obtain the loffe-Larkin composi- 
tion rule which states that the inverse of the superfluid 
density ps — ns/m is given by adding the inverses of the 



fermion and boson contributions: 



(13) 



where pf « xt and pf « J(l - aT) with a « A^^ ^ 
Note that the linear temperature dependence comes from 
thermal excitations of the nodal fermions. Expanding 
Eq. ^|at small T, we obtain 



Ps{T) « pM 



pfio) 



aT. 



(14) 



On the other hand, the quasiparticle dispersion is given 
by that of the fermions and we can identify vf/va in Eq. 
^as being proportional to a J ^ apf{0). Comparison of 
Eq. 1141 with Eq. ^ results in the simple expression 



pM_ 

pf (0) ' 



(15) 



in particular a oc x{t/J) for xt < J. [23| We did not 
keep track of the numerical coefficients. If we assume 
that at full doping {x = 1) a should approach one (the 
BCS value), Eq. lfTC|l suggests the form 



a{x) 



:(1 



(16) 



which for x* — 0.09 produces a qualitatively good fit 
of our results for a (see Fig. I^J. Note that the order- 
of-magnitude estimates above gives x* J/t, and our 
best-fit value of x* is several times smaller. 

The shape and the magnitude of the current renor- 
malization parameter a as a function of doping is cur- 
rently a topic of much experimental work. While there 
still remain large uncertainties in the experimental data, 
two useful comparisons can be made to our work. In 
Ref. 0, the measurements of thermal conductivity and of 
the penetration depth were used to determine the ratio 
vf/va and the superfluid density. Combining those data 
resulted in the the values of a = 0.66 and a = 0.68 for 
optimally doped samples of BSCCO and YBCO, respec- 
tively. Those numbers are in close agreement with our 
results around x — 0.15 doping. 

Our results also qualitatively agree with the decrease 
of a as the doping decreases in underdoped YBCO, as 
reported in Refs. 0, 0. However, earlier data on YBCO 
films indicated that the linear T slope of ng/m is rela- 
tively insensitive to dop ing over a broad range of critical 
temperatures. [2^ l23| Since vf/va decreases with de- 
creasing X ll^l, this trend is in disagreement with Fig. [31 
and Eq. (|16l) . It will be desirable to have single crystal 
data for (T) / m over a broad range of x to settle this 
point. 

On the theoretical side, the recent study of a U{1) 
slave-boson theory with spinon-holon binding [j^ l pre- 
dicted a sub-linear dependence of a on the doping. The 
sublinear form of a{x) disagrees with our proposal H16|l . 
but is also consistent with our numerical results shown 
in Fig. [3 

We further look at a along the Fermi surface for a 
given doping. We are interested in how the shape and 
size of this curve changes as a function of doping. While 
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FIG. 6: Plot of |j| and a = |j| /«f, upper and lower plots re- 
spectively, for the t-J model along the Fermi surface for dop- 
ings of 0.01,0.05,0.09,0.17 and 0.25 (increasing magnitude). 
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we expect the integral of the whole curve to decrease 
with X, there are several possibilities for how this could 
occur. Two such scenarios are that the curve decreases in 
magnitude everywhere along the Fermi surface uniformly 
as doping decreases or that a is small almost everywhere 
along the Fermi surface but that there is a region of large 
a around the nodal point whose width increases with 
X. Our numerical results indicate the former of those 
scenarios. In the upper graph of Fig. we plot the 
t-J model current magnitude along the Fermi surface. 
Due to the low finite resolution of our system, we just 
calculate the current at those points nearest to the line 
connecting (7r,0) and (0, tt). To calculate a, we use the 
value of vp found by Yunoki et al. 3J along the nodal 
direction for a given doping and accordingly renormalize 
the tight binding dispersion to calculate Vk along these 
points. We plot a obtained in this way in the lower plot 
of Fig. El We see that a is approximately flat along the 
Fermi surface. 



III. QUASIPARTICLE WEIGHT AND 
OCCUPATION NUMBER 



In this section we examine the quasiparticle weight Z 
and the occupation number n as a function of momentum 
and doping. The quasiparticle weight is defined as in 



FIG. 7: Contour plots of .Z+(k) in the t-J model for dopings 
0.03 (upper plot) and 0.21 (lower plot), 10 and 68 holes in 
a 18x18 system, respectively. The thick dashed lines denote 
the Fermi surface of the unprojected wave function, the big 
solid dot marks the position of the node. Crosses indicate 
data points used in plotting (available values of the k vector) . 
In the upper plot, the contour lines correspond to Z+ = 0.01, 
0.02, . . ., 0.09 (left to right - the maximal value of Z+ is 0.09). 
In the lower plot, the contour lines are Z+ = 0.05, 0.10, . . ., 
0.35 (with the maximal value Z+ = 0.38). 



Fermi liquid theory and gives a measure of how close our 
trial wave function quasiparticle is to being a free electron 
(or a free hole). 
We begin with the definitions. 



^+(k,a) 



Z_(k,(7) 



;^Ex|<J4-(3s) 



(*EX|*EX)(*GS|*GS) ' 
(^Ex|ck,^|^Gs)|^ 

(«'ex|«'ex)(^'gs|^'gs) ' 



(17) 

(18) 



where |4'ex) also carries momentum k, and the electron 
operators are normalized as {ck.o-, Ck ^.j = 1. In the BCS 
theory (without Gutzwiller projection), = and 
Z_ = v^. Along the nodal diagonal (where the gap van- 
ishes). 



1 outside the Fermi surface and inside and 



TJk is the opposite. If one defines Z = Z+ + Z-, then 



m 
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the BCS model Z = \ everywhere. 

In the projected wave functions, these simple expres- 
sions do not apply. However, some of the properties of 
the spectral weights Z+ and Z_ may be still proven. Us- 
ing the identity 

Pdc\Pd = Pdc\ (19) 

and the d-wave symmetry of the gap, we can prove that 
along the nodal diagonal, 

Z+ = on (0,0) - (fcF,fcF), (20) 
Z_=0 on (fcF,fcF) - (7r,7r), (21) 

where (/cf, ^f) is the nodal point. Furthermore, Z+ may 
be rewritten as the ground-state expectation value 

Z+(k,a) = (ck,,F,5C^)GS. (22) 

We note that 

(q.<t4^)gs = (c,;,<T^'c],f^)GS (23) 
= (c,^<,f4^)gs 1^3- (24) 

Thus (k, cr) is further related to the occupation num- 
ber 

n\i,a = (Ck^^Ck.o-) (25) 

as 

Z+(k,fT) = l±^-nk,.. (26) 

This relation has also been given by Yunoki in Ref. 
In particular, from this relation follows the upper bound 
on the spectral weight Z^: 



and, as a consequence, the same upper bound applies to 
the nodal spectral weight Z„odai = ■^i_(^£,+ ^i + e) 
studied by Paramekanti et al. in Refs.l2.llOl 

The quasiparticle weight Z^ requires a more compli- 
cated Monte Carlo calculation, since it cannot be rewrit- 
ten as a simple ground-state average like (f^ . We there- 
fore restrict ourselves to discussing only the quasiparticle 
weight Z+ in this paper. 

The above relations for Z+ have been derived for the 
fully projected wave function. If we perform a rotation 
to the Hubbard model, the relations (|22() and 1261) no 
longer hold, and the upper bound (|27f) cannot be proven. 
Specifically, for the Hubbard model we keep the same 
definitions of the quasiparticle weights and the occupa- 
tion number (|17l) . (|18l) . (|25|l . but with the ground and 
excited states rotated from the fully-projected state to 
the Hubbard- model state by the unitary rotation e"*"^, 
as explained in the previous section. Then the lowest- 
order Hubbard-model corrections to Z^ and to may 
be easily computed as 

Z^ = Z++2Re (ck,.PD[z5, 4_J) (28) 




0.2 0.4 0.6 0.8 1 

kx, units of TT 




0.2 0.4 0.6 0.8 1 

kx, units of TT 




0.2 0.4 0.6 0.8 1 

kx, units of TT 



FIG. 8: Upper plot: contour plot of ^+(k) (or, equiv- 
alently, of n(k)) in the t-J model at the doping 0.12 (24 
holes in the 14x14 system). The isolines correspond to 
Z+ = 0.03, 0.06, . . . , 0.27 [n(k) = 0.53, 0.50, . . . , 0.29 respec- 
tively, left to right] . The maximal value of Z+ (near the node) 
is 0.27. Middle plot: Z^ in the same system rotated to the 
Hubbard model (to the first order in t/U = 1/12), Eq. ij^ . 
The isolines correspond to Z+ = 0.05, 0.10, 0.35. The 
maximal value of Z^ (near the node) is 0.37. Bottom plot: 

in the same system rotated to the Hubbard model (to 
the first order in t/U = 1/12), Eq. The isolines are 

= 0.8, 0.7, . . . , 0.2 (left to right). The Fermi surface, the 
nodal point and the positions of data points are denoted in 
the same way as in Fig.Q 
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0.4 



0.3- 



0.1 - 



0' 




(0,0) (Ji/2,ji/2) (Ji,ji) 

FIG. 9: Plot of and , the solid and dashed lines re- 
spectively, for an 18x18 system with 40 holes, x = 0.12. The 
data are calculated on the discrete grid points closest to the 
diagonal. 

and 

= nk - 2 Re (ck,.[i^, c[ J) . (29) 

Even though the two expressions look nearly identical, 
the correction to rik does not contain an intermediate 
projector Pq and, because of that, has a very different 
structure than that to Z+ . The relation between Z+ and 
TT-k H26|) no longer holds for the Hubbard model, as we 
shall see below. 

Figs. El and m show Z+(k) for x = 0.03, 0.12, and 
0.21. As noted in Eq. (|20|l . Z+ is zero along the diag- 
onal between (0,0) and (kp^^F) and jumps to a finite 
value Znodai at the nodal point. The projected wave 
function inherits the essential singularity of at the 
nodal point from the underlying BCS wave function. The 
value of 2'nodai has been studied as a function of doping 
by Paramekanti et al. in Refs. 0, The doping de- 
pendence of .^nodai (Fig. 2 of Ref. and Fig. 6 of Ref. 
llOr) is qualitatively similar to that we find for a (Fig.[SJ): 
both 2'nodai and a decrease to zero with decreasing dop- 
ing, with a strong upward curvature. However we are not 
aware of any a priori relation between a and Z^odai- 

Using iCni), the plots of Z+(k) may also be interpreted 
as those of rik (see, e.g., the upper plot in Fig.|SJ|. Note 
the region of depression in rik just outside the Fermi sur- 
face, which resembles a hole "pocket." The existence 
of this pocket may already be inferred from the non- 
monotonous behavior of rik along the zone diagonal found 
in Refs. but the full k dependence shown in Figs.Q 
and|Slgives a more complete picture. Remarkably, a sim- 
ilar "pocket" structure has been found in the U (1) slave- 
boson model with spinon-holon-binding by Ng |25| . The 
pocket is more pronounced at lower dopings and appears 



consistent with the bending of the "effective Fermi sur- 
face" defined in Section ^ from nj , Eq. However, to 
define a meaningful Fermi surface from the quasiparticle 
spectral weight, one needs an access also to the spectral 
weight Z- which goes beyond the scope of the present 
paper. 

If we include the t/U correction from the rotation to 
the Hubbard model, the pocket structure in nff disap- 
pears, see our Fig. |S1 (lower plot) and Refs. 0, On 
the other hand, including the t/U correction to Z+(k) 
preserves the pocket structure, see Fig. |S1 (middle plot). 
This shows that in the Hubbard model the relation l|2t)|) 
between Zj^ jk) and no longer holds. 

In Refs. a, M it was reported that the rotation to the 
Hubbard model does not change the magnitude of the the 
jump in rik at the Fermi surface. Our results on rik and 
confirm this statement, however the spectral weight 
defined as Z+ increases when rotated to the Hubbard 
model. In Fig. ^ we show Z+ and Z^^ along the nodal 
diagonal for an x ~ 0.12 system. 



IV. CONCLUSION 

In this paper, we have analyzed the properties of the 
excited states in the t~J and Hubbard models in the 
framework of Gutzwiller-projected variational wave func- 
tions. The quantities of main interest are the renormal- 
ization of the current and spectral weight of the quasi- 
particles. Both those renormalizations decrease with de- 
creasing doping and exhibit strongly non-BCS behavior. 

The renormalization of the quasiparticle current allows 
us to define the effective Fermi surface as a crossover 
region between the electron- and hole-supported current. 
We observe that such a Fermi surface bends outwards 
in the (0, tt) regions - an effect normally ascribed to the 
t' hopping term. At the same time, the total current is 
renormalized approximately uniformly along the Fermi 
surface. 

The renormalization of the quasiparticle spectral 
weight, on the other hand, is peaked near the nodal point 
and exhibits a pocket-like structure. This pocket feature 
is more pronounced at lower doping values. 

Comparing the results for the t-J model and for the 
Hubbard model (to the lowest order in the t/U correc- 
tion), we find that the rotation to the Hubbard model 
does not qualitatively affect the renormalizations of the 
quasiparticle spectral weight and of the current. 

We would like to thank Professor T.K. Lee and Pro- 
fessor Castellani for their discussions and help on this re- 
search. C.P.N, and P.A.L. acknowledge support by NSF 
grant DMR-0517222. 
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